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Abstract — This paper introduces a space-continuous force- 
based model for simulating pedestrian dynamics. The main 
' interest of this work is the quantitative description of pedestrian 
movement through a bottleneck. Measurements of flow and 
density will be presented and compared with empirical data. 
The results of the proposed model show a good agreement with 
empirical data. Furthermore, we emphasize the importance of 
volume exclusion in force-based models. 

I. Introduction 

Recently pedestrian dynamics has been gaining increas- 
ing interest. One focus of this research area is the security 
of people, e.g., trying to optimize evacuation processes by 
minimizing evacuation time. To achieve this "benchmark" 
architects and civil engineers need to have an idea about 
the minimal necessary width of exit doors, their optimal 
placement, length of escape routes, etc. Usually legal reg- 
ulations and other descriptive specifications in handbooks 
provide requirements, which are in general not flexible enough 
] for complex buildings. However, simulations of pedestrian 
I dynamics offer the possibility to analyze and understand 
the dynamics of pedestrian streams in the building to be 
designed. This is useful to ease decisions e.g., dimensioning of 
emergency doors, and make them more realistic and adapted 
to different architectures. 
] In this work, we address the possibility of describing 
! quantitatively the pedestrian dynamics, by proposing a space- 
continuous model. In section 2 the theory of mathematical 
modeling of pedestrian dynamics is briefly discussed. In 
section 3 our model is introduced. Then in section 4 numerical 
results of pedestrian flow through bottleneck and density mea- 
surements inside and in front of the entrance to a bottleneck 
will be presented. Section 5 gives concluding remarks. 

11. Mathematical models 

There is a wide range of mathematical models to simulate 
pedestrian dynamics. Generally, these models are subdivided 
into macroscopic and microscopic models Q. In macroscopic 
models the system is described by variables like crowd-density 



and flows of continua, whereas microscopic models consider 
the movement of individual persons separately. Microscopic 
models can be subdivided into models which are discrete or 
continuous in space. Widely-used discrete models for pedes- 
trian dynamics are Cellular Automata models (CA), which 
describe phenomena in space-time by assigning discrete states 
to a grid of space-cells. These cells can be "occupied" by a 
pedestrian or "empty". Thus the movement of pedestrians in 
space is done by passing them from cell to cell in discrete 
time by simple rules that reflect e.g. psychological aspects of 
the motion. 

Space-continuous force-based models determine the contin- 
uous movement of well distinguished individuals from their 
desires to reach certain targets and from influences of the space 
geometry as well as actual states and positions of the other 
individuals. Inspired by the field theory of Lewin ||4l, force- 
based models assume that not only physical forces do affect 
the movement of pedestrians but similarly do social forces. A 
mathematical description of such social forces that determine 
the movement of pedestrians was first given by Helbing and 
Molnai- in 1995 12J. 

Though force-based models pose quite some computational 
problems that do not occur in CA, their investigation appears 
to be worthwhile, since they permit higher resolution of 
geometry and time. 

Force-based models take Newton's second law of dynamics 
as guiding principle. The movement of each pedestrian is thus 
described by 



F^=T. + E + Ft'' = m,^. (1) 
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Here F['^'^ denotes the repulsive force exerted by pedestrian 

> 

j on pedestrian i. F[^^ is the repulsive force emerging from 

borders (e.g. walls) and Ff^"" is a driving force, nii is the 
mass of pedestrian i and N is the number of pedestrians 
acting on pedestrian i. Pedestrians try to avoid collisions and 




Fig. 1. Direction of the repulsive force acting on a pedestrian 



on the inverse of their distance. Additionally, relative velocity 
between two pedestrians is taken into account and assumed 
to influence the strength of the repulsive force. Moreover, 
we assume proportionaUty between the repulsive forces acting 
on a pedestrian i and its desired velocity. Thus a pedestrian 
tending to move with high velocity "feels" intense repulsive 
actions. 

The repulsive force, as a dimensionless relationship between 
the above mentioned quantities, is defined by 



contact with other pedestrians and objects by changing their 
direction. The repulsive forces are introduced to model this 
avoidance of physical contact. They can be social forces |2|, 
m, IIT2I as well as physical forces HI, |j5J, where physical 
forces come into play mainly to prevent overlapping with other 
pedestrians and with obstacles. Several different forms for 
these repulsive forces have been proposed. The driving force 
pdrv represents the intention of a pedestrian to move towards 
a given destination. 

The set of equations ([T]i for all pedestrians results in a 
high-dimensional system of second order ordinary differential 
equations. The time evolution of the positions and velocities 
of all pedestrians is then obtained by numerical integration. 

Most force-based models describe the movement of pedes- 
trians qualitatively well; Self-organisation phenomena e.g., 
lane formations Q, |[T|, lfT2l . oscillations at bottlenecks ||2l, 
im, "faster-is-slower" effect f3l, fSl, and clogging at exit 
doors im, lfT2l are reproduced. These achievements indicate 
that the models could rehably describe pedestrian dynamics. 
However, a qualitative description is not sufficient if one 
has to make reliable statements about critical processes, e.g., 
emergency egress. 

Moreover, rather often software implementations of the model 
do not rely on one sole approach. Especially in high density 
situations the simple numerical treatment has to be supple- 
mented by additional techniques to obtain sensible results. Ex- 
amples are restrictions on state variables and sometimes even 
totally different procedures which replace the above equations 
of motion ([T]i to avoid overlappings among pedestrians 131, 
lfT2j and negative and high velocities [2J. 

Necessarily a model has to describe pedestrian dynamics in 
a quantitative manner, e.g., by reproducing the fundamental 
diagram, the flow through bottlenecks and the density inside 
and in front of the entrance of a bottleneck. Additionally, for 
further investigation, it is advantageous that it uses a single 
set of forces. 

In the next section, we propose such a model, which is 
based on a single set of formulae. Furthermore the model 
incorporates free parameters which allow adaption of the 
model to quantitative data. 

III. A Modified Centrifugal Force Model 

Pedestrians are modelled as circles moving in two- 
dimensional space, see Figure [T] We model the repulsive 
force exerted by pedestrian j on pedestrian i as depending 
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where 
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with R,j = Rj - R„ (3) 



and 



and 



V^J = - Vj) -^j + m- Vj) ■ ^|], (4) 



Vi ■ eij+ I Vi ■ ejj 
II ^11 



(5) 



Here Vf' denotes the desired velocity of pedestrian i and rrii 
its mas^. Vij is the projection of the relative velocity of 
pedestrian j and pedestrian i onto the direction of Cij, \\ Rij \\ 
is the distance between i and j, Kij is a coefficient that reduces 
the action-field of the repulsive force to the angle of vision 
of each pedestrian (180°), and is the normalised direction- 
vector between i and j. \\ . \\ denotes the Euclidian norm in 
and |.| the absolute value in JR. 
By means of the parameter i' the strength of the force can 
be adjusted while the diameter of pedestrians Di{\\ Vi ||) 
depends linearly on the velocity. This incorporates the dynamic 
space requirement of pedestrians, modeling the fact that faster 
pedestrians require more space than slower pedestrians, due 
to increasing step sizes 



A = Da + Db II V, II 

with free parameters Da and Df,. 

We use the driving force as defined in 
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(6) 
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with Vi the velocity of pedestrian i, e" its desired direction 
and r a time constant. 

The above model has been developed from the Centrifugal 
Force Model (CFM) of Yu et al. U_2 |. Their expression for the 
repulsive force reads as follows: 



-miKij- 



R^ 



(8) 



In all simulations the mass is set to one for all pedestrians. 



with the quantities Kij, Vij, Rij and Eij as defined in 
Equations Q, and Q. 

When pedestrians come near to one another the repulsive 
force is expected to grow. In the CFM, this is in general not 
the case due to the relative velocity term in the nominator 
of the repulsive force (HJ. If both the distances between two 
pedestrians and their relative velocities are small, their quotient 
does not grow as expected. Consequently the repulsive force 
does not become high enough to prohibit overlappings. Intro- 
ducing the intended speed in the nominator of the repulsive 
force (O we eliminate this side-effect. Furthermore, there is 
a compensation between repulsive and driving forces at low 
velocities, which damps oscillations. Due to these changes we 
can do in most cases without the extra Collision Detection 
Technique (CDT) which takes over control in 1,12] in case of 
formation of dense crowds. 

IV. Numerical results 

The solution of the initial value problem (HJ was done in 
all simulations using a fourth-order Runge-Kutta scheme with 
a fixed-step size of At = 0.01 s. The repulsive interaction be- 
tween two pedestrians is neglected when the distance between 
them is larger than two meters. 

In order to investigate the influence of our modifications, 
different model approaches were simulated. Then the results 
were compared with empirical data of pedestrian flow through 
a bottleneck. 

The flow of 60 pedestrians through the bottleneck situation 
as described in |8| was simulated. The width of the bottleneck 
was changed from 0.8 m to 1.2 m in steps of 0.1 m. The free 
parameter v in Equation (|2]i is set to 0.28 m/s for the repulsive 
force among pedestrians and to 0.4 m/s for the repulsive 
force emerging from obstacles. For the free parameters in 
Equation ^ we set 

Da = 0.2 m and Db = 0.2 s. 

Several parameter values were tested. With this parameter 
set the results of the simulations are in good agreement with 
the empirical flow measurements presented in |7|, see Fig- 
ure [2(a)] For comparison the flow through the same bottleneck 
was measured by the original CFM without CDT. The results 
presented in Figure |2(b)| show very high values of the flow. 
These are possible, because CDT does no more care for 
volume exclusion. 

Actually, control of bottleneck-flow seems to be dominated 
by the CDT in the CFM-approach. To back this hypothesis we 
measured pedestrian flow through the bottleneck with a CFM 
like model without repulsive forces, managing collisions with 
CTD, only. The gained values are in the range of experimental 
data, see Figure |2(c)| Volume exclusion of pedestrians seems 
thus to be a favorable feature of successful models. 

A second validation of our modifications of CFM comes 
from measurements of density inside the bottleneck as well 
as in front of the entrance to the bottleneck, see Figure |3] 
The density in front of the entrance to the bottleneck is 
presented in Figure |4(a)| The results are in good agreement 
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(a) Flow measurement with the modified CFM without CDT 
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(b) Flow measurement with original CFM without CDT. Without CDT the 
flow values are unrealisticly high. 
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(c) Flow measurement with CDT without repulsive forces. 



Fig. 2. Flow measurement with tree different models in comparison with 
empirical data (7). 



with the experimental data in ifTOl . Additionally, the measured 
density values inside the bottleneck are in accordance with the 
published empirical results in |8|, see Figure [4(b)] 
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Fig. 3. Screenshot of a simulation: The density of pedestrians is measured 
by counting tiie number of pedestrians that are located inside one of the 
measurement rectangles at each integration time. The area of the appropriate 
rectangle divides the result. 
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(a) Density in front of the entrance to the bottleneck 
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(b) Density inside the bottleneck 



oscillations and overlappings can still occur Determination of 
uniformly usable parameters is a difficult endeavour; cf. the 
parameter estimation in SimWalk 

V. Conclusion 

The proposed space-continuous force-based model quantita- 
tively describes the movement of pedestrians in 2D-space in an 
almost satisfactory way. Besides being a remedy for numerical 
instabilities in CFM the model simplifies the approach of Yu 
et al. since we can dispense with their extra CDT without 
deteriorating performance. The implementation of the model 
is straightforward and does not use any restrictions on the 
velocity as was the case e.g., in 111. 

Simulation results show good agreement with empirical 
data. The model contains natural free parameters which can 
further be tuned to adapt the model to certain scenarios. 
Introducing density-dependence of the strength of the repulsive 
force promises further ameliorations of the model. 

Furthermore we have provided some evidence that volume 
exclusion is an important issue in describing pedestrian dy- 
namics. 
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Fig. 4. Density measurements: The simulation results (blue lines) are in 
good agreement with the empirical data presented resp. in 1101 and (8). 

Though the model is already a useful tool to describe 
quantitatively the bottleneck-situations as defined in [81, it is 
not yet fully satisfactory. Depending on the specific situation 



